A1 = readmatrix('-10V_alongQ.txt')


x1 = A1(1:90,1)/360*2*pi
y1 = A1(1:90,2)*1e12-0.06503*x1

a = -0.1362
b = -0.07457
c = -0.1085
d = -1.479

f1 = c*sin(2*x1)

hold on
plot(x1,f1,'Linewidth',2,'Color',[0.8 0 0])
scatter(x1,y1,10,[0 0 0.8],'filled')
xlabel('Polarization')
ylabel('Photocurrent [pA]')
xlim([0, pi])
ylim([-0.15,0.15])
xticks(0:pi/4:pi)
yticks(-0.1:0.1:0.1)
xticklabels({})
yticklabels({})
%set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 200, 80])